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Abstract. We find a normal form which describes the high energy dynamics 
of a class of piecewise smooth Fermi-Ulam ping pong models; depending on 
the value of a single real parameter, the dynamics can be either hyperbolic 
or elliptic. In the first case we prove that the set of orbits undergoing Fermi 
acceleration has zero measure but full Hausdorff dimension. We also show that 
for almost every orbit the energy eventually falls below a fixed threshold. In 
the second case we prove that, generically, we have stable periodic orbits for 
arbitrarily high energies, and that the set of Fermi accelerating orbits may 
have infinite measure. 



1. History and introduction 

In this paper we study the dynamics of piecewise smooth Fcrmi-Ulam ping pongs; 
Fermi and Ulam introduced such systems as a simple mechanical toy model to 
explain the occurrence of highly energetic particles coming from outer space and 
detected on Earth (the so-called cosmic rays, see [TSIIIS])- The model describes the 
motion of a ball bouncing elastically between a wall that oscillates periodically and a 
fixed wall, both of them having infinite mass. Fermi and Ulam performed numerical 
simulations for the model and consequently conjectured (see |25j ) the existence of 
orbits undergoing what is now called Fermi acceleration, i.e. orbits whose energy 
grows to infinity with time; we refer to such orbits as escaping orbits. Several years 
later, KAM theory allowed to prove that the conjecture is indeed false. Namely, 
provided that the wall motion is sufficiently smooth, there are no escaping orbits 
because invariant curves prevent diffusion of orbits to high energy (see POIIMII^ ). 
It was not many years (see [2H|) before the existence of escaping orbits was proved in 
some examples of piecewise-smooth motions; it is worth noting that these examples 
were essentially the same that Fermi and Ulam were forced to investigate in their 
numerical simulations, due to the relatively limited computational power they could 
us^ In this paper we study a more general class of piecewise smooth motions and 
we investigate existence and abundance of escaping orbits in this setting. 

Our main result is that, for all possible wall motions having one discontinuity, 
there is a parameter A which allows to describe the dynamics of the ping pong for 
large energies. Moreover, there exists a sharp transition so that for A S (0, 4) the 
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system looks regular for large energies while for A ^ [0, 4] the system is chaotic 
for large energies (see Figure [2] in Section [2]) . Similar phenomena happen in a 
wide class of piecewise smooth mechanical systems which for large energies can be 
viewed as small (non-smooth) perturbations of integrable systems such as, for 
example, the impact oscillator |15j . However, in order to demonstrate the methods 
and techniques in the simplest possible setting, we restrict our attention to the 
classical Fermi-Ulam model. 



2. Results. 

We consider the following one-dimensional system: a unit point mass moves 
horizontally between two infinite mass walls, between collisions the motion is free, 
so that kinetic energy is conserved, collisions between the particle and the walls are 
elastic. The left wall moves periodically, while the right one is fixed. The distance 
between the two walls at time t is denoted by £(t) which we assume to be strictly 
positive, Lipshitz continuous and periodic of period 1. It is convenient to study the 
orbit only at the moments of collisions with the moving wall. Let t denote a time of 
a collision of the ball with the moving wall, since i is periodic we take i e T = M/Z. 
Let G M be the velocity of the ball immediately after the collision. Introduce the 
notation A = T x M; the collision space is then given by 

M = {(t,w) e A s.t. V > ~l{t)}. 

We can thus define the collision map f : A4 ^ M 

(1) fitn,Vn) = {tn + St(tn,Vn),Vn - 2l{tn + dt{tn,Vn))) ^ (t«+l,W„+l) 

where, for larg^u, the function St solves the functional equation 

(2) Stit,.)^'-^^ttl(l±m^, 

V 

It is a simple computation to check that the map / preserves the volume form 
LJ = {v + £{t))dt A dv. Throughout this work we assume £ to be piecewise smooth 
with a jump discontinuity at t = only. Define the singularity line 5 C as 
S — {t = 0} and R G Ai as the infinite strip of width O {v~^) bounded by S and 
fS; introduce also R — f^^R- 

As a first step to study the dynamics of the mapping / we describe the first 
return map of / to the region R, which will be denoted hy F : R R. Our 
main result is a normal form for F for large values of v: introduce the notation 
£o = £{0), £^ = i(0±) and similarly for aU derivatives. Define A = J£o{£+ - t) 
and Ai = \j'^£l(£+ - 't) where J is given by 

(3) I £-^{s)As. 



We introduce a useful shorthand notation. Let S C*(v4 C A); then we use the 
notation ip — Og to indicate that v'^ip is bounded for sufficiently large v and 

the same is true for all derivatives up to order s included. For our analysis it is 
important to ensure that all sub-leading terms vanish sufficiently fast for v ^ oo 
along with all partial derivatives up to the fifth order. 



^ large here means that the ball bounces off the fixed wall before the next collision with the 
moving wall 
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Theorem 1. There exist smooth coordinates (t,I) on R, such that the first return 
map of f on R is given by 

F{t, I) = F{t, I) + i^i (r, /) + r (r, /) 

where F{t, I) ~ (f, /) with 

f = T-I modi, / = / + A(f-l/2), 

Fi is a correction of order O {l^^) of the form 

^i(r,/)=/-i(0,Ai((f- 1/2)2 -1/12)) 

and r = O5 {l^^) ■ Finally u = dr A dl. 

Consequently, up to higher order terms, F coincides with F, where F is Z^- 
periodic in appropriate "action-angle" variables and moreover dF = A is constant. 
Thus F covers a map F : — )■ T^; the map F is known in the literature as the 
"sawtooth map" or the "piecewise linear standard map" and it has been the subject 
of a number of studies, see e.g. [21 131 SI El 13 [57]. Notice that we have 

Ti-{A) = 2 - A. 

Accordingly, dF is elliptic if A € (0,4) and it is hyperbolic otherwise. 
Example 2.1. Consider the case where velocity is piecewise linear, that is 
(4) £A^B{t) = B + A{{t nlodl)-l/2)^ 

This is one of the cases which have been numerically investigated in [25^. We 
can choose the length unit so that B = 1. In this case £ is positive for all t for 
A > —4. Remarkably, we can obtain an explicit expression for A{A), that is, for 
J {A). Namely 

JiA) = / (1 + A{s - l/2f)-^ds = 21^1-1/2 / (1 + sgnA • a^y^da 

Jo Jo 

where a = \A\^/^{s — 1/2). Performing the integration we obtain 

-..^^ 2 ^ /(|A|-VV2)logg^ if - 4<A<0 
A + A \|A|-i/2arctan(|A|i/2/2) if A > 0. 

Recall that, by definition, A{A) ^ -2A{1 + A/A)J{A). In particular, wc find that 
A(-4) = 4 and 

1 2 + lAI^/^ 
lim A'(^) = lim -log ' ' ^ ^ +00. 

The graph of dependence of A on ^ is shown on Figure [Tj It shows that the 
dynamics is hyperbolic for A £ (—4, —a) where a ~ —2.77927 and A > and it is 
elliptic for remaining parameter values. 

Below we discuss the implications of the dichotomy between hyperbolic and 
elliptic regimes. Using the results of [10' we obtain the following 

Theorem 2. //|Tr(y4)| > 2 then F is ergodic, mixing and enjoys exponential decay 
of correlations for Holder observables. 
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Figure 1. Graph of A as a function of A, where Ia.b is given by 
Q with B = 1; the shaded area denotes the eUiptic regime A G 
(0,4). 

On the other hand if |Tr(A)| < 2 then F is not ergodic. Namely, in this case, 
_F is a piecewise isometry for the appropriate metric. Hence if p is a periodic point 
of F, then a smaU baU around p is invariant by the dynamics. See Figure [2] for an 
example of phase portrait of F in the two cases. 

Note that if p is periodic with period N for F, it need not be periodic for F. In 
fact we have 

F^p - p + (0, n) 

for some n S Z. If n > we say that p is a stable accelerating orbit; if n < we say 
that p is a stable decelerating orbit; finally if n = then p is periodic for F. 

Consider for example the case = 1: we can find a periodic orbit (|,0); if 
furthermore A > 2, we have a stable accelerating orbit (0, 5 + ^) and stable 
decelerating orbit {0,^ — 

To analyze periodic points we can use the duality between the accelerating and 
decelerating periodic orbits. We have F = Ta o G where 

(5) G{t,I) = {t-I modi,/), rA(r,/) = (r,/ + A(r-l/2)). 

On the other hand F^^{f, I) — (r, /) with 

/ = /- A(f- 1/2), T = f + I. 

Introducing ct = 1 — t wc rewrite the last equation as 

/ = / + A(cr-l/2), cr = CT-/. 

In other words if J denotes the involution J(t, /) = (1 — r, /) then 

(Ta o G)"1 = J oGoTao J = 

= o G o Ta o J = (Ta o J)-^ o (Ta o G) o (Ta o J). 

The existence of periodic orbits for other small periods is summarized in table [l] 
(here we use parameter 9 such that Tr(dT) ~ 2 cos 6*, that is, dF is conjugated to 
a rotation by 9) . 
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Figure 2. On the top: phase portrait of a single orbit of the map 
for A = —0.3. On the bottom: phase portrait of selected orbits 
of the map F for A = 0.32. Notice the strong prevalence of elliptic 
behavior; the "chaotic" region is given by forward and backward 
images of the singularity line. 
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N periodic 



accelerating/ decelerating 



1 (0,7r) 

2 (0,7r/2) 

3 (0,V3) 

4 (0,7r/4) 

5 (0,^5) 



(3/47r,^) 
(7r/2,3/5^) 

(5/67r,^) 
(3/7^,7r/2) 

(7/87r,7r) 



(7r/2,7r) 



6 (0,7r/6) U (7r/2,7r) 



7 (0,^/7) 

8 (0,7r/8) 



Table 1. Numerically observed ranges of parameters for which 
there exist a periodic or accelerating/decelerating orbit of "pe- 
riod" N; ioi N > 9 we could find intervals bounded by irrational 
multiples of tt 



Remark 2.2. We believe that stable escaping orbits should exist for arbitrarily small 
positive values of A, i.e. for each Aq > there exists a < A < Aq such that 
the map F admits stable escaping orbits. However their "period" will necessarily 
to grow to infinity as A — >• 0+; the smallest value of A for which we were able to 
find a stable escaping orbit is A = 0.0916346, for which we numerically obtained a 
period 501 stable escaping orbit. 

We return now to the problem of energy growth; introduce the set of escaping 
orbits 



Theorem 3. If\Tr{A)\ < 2 then 

(a) there exists a constant C so that for each v sufficiently large there exists an 
initial condition (to,Vo) such that 



If additionally Ai ^ 0, then the same result holds for an adequately small ball 

around the point (to.wo)- 
(b) If F has a, stable accelerating orbit then mes{£) = oo. 

In particular mes (f) =oo if Ae (^3,4). 

Theorem 4. //|Tr(^)| > 2 then 

(a) mes{£) = 0. 

(b) there exists a constant C such that almost every orbit enters the region v < C. 
Moreover denote by T the first time velocity falls below C. If we fix the initial 
velocity 1 and let the initial phase be random then ^ converges to a stable 

random variable of index 1/2, that is, there exists a constant D such that 



The proof of second part of the last theorem relies on the following result which 
is of independent interest. 

Theorem 5. Fix the initial velocity wq 1 and let the initial phase be random 
then Fix < a < 1 < 6. Consider the process defined by B^o(i) = if vlt is 



£ = {{tQ,vo) : Vn oo}. 



(6) 



C < Vn < Cv for all n e Z. 
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an integer with linear interpolation in between which is stopped when velocity goes 
above bv^ or below av^. Then, as Vq oo, B''"(i) converges to a Brownian Motion 
started from 1 and stopped when it reaches either a or b. 



Remark 2.3. Note that 'B'"''{t) is equal to ' only if v^t is an integer. It seems 



more natural to use this formula for all values of t however this would lead to a 
different limit since as we shall see in the next section the ratio v{n + ^)/v{n) has 
oscillations of order 1 while v{n + l)/v{n) — 1 is of order l/v{n). 

Theorem [5] makes Theorem |4] plausible since the time the Brownian Motion 
reaches a certain level has a stable distribution of index 1/2. However some work is 
needed to deduce Theorem [4] from Theorem [5] since the proof of Theorem [5] relies on 
a perturbative argument near v = oo whereas Theorem [4] requires to handle small 
velocities as well since w(T) < C. 

Theorem |4] shows that the set of escaping orbits has zero measure so it is natural 
to ask about its Hausdorff dimension. The next result extends the work [T3] where 
a similar statement is proven for a smooth model of Fermi acceleration. 

Theorem 6. //Tr(A) > 2 then HD{£) = 2. 

In other words, even though the set of escaping points is small from the measure 
theoretical point of view, it is large from the point of view of dimension. 

3. The first return map 

If £ were a smooth function, KAM theory would allow to conjugate the dynam- 
ics of / for most initial conditions with large values of v to the dynamics of the 
completely integrable map g : T x M+ — > T x M+ 

g:{^,J)^{^ + J-\j). 

Consider the vertical line 5' C T x K+ given by S' = {d = 0}; let moreover R' be 
the infinite strip of width O {J~^) bounded between S" and gS' i.e. 

i?' = {0 < z9 < J-^}. 

As a preparatory step we study the first return map of g to the region R' . 

Proposition 3.1. Let t — and consider coordinates {t,J) on R' . Then the 
first return map of g to the region R! is given by the map G defined in 

Proof. Let k = [ JJ and J = k + J. We claim that 

'g^{T,J) ifT<l~J 



G(t, J) 



g^ {t,J) otherwise. 



In fact, we can check by simple inspection that, denoting g''{{},J) = J) we 
have 

k + J J J 

which implies our claim. □ 

In our systems, £ is only piecewise smooth, consequently we expect to be able to 
define action-angle coordinates outside R only. 
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Lemma 3.2 (Approximate reference coordinates). There exists a smooth coordi- 
nates change h : (t, v) i-> (-(9, J) such that if {t, v) ^ R, h conjugates the collision 
map f to the reference map g up to high order terms 

(7) g ~ ho f o h"'^ = {r^,rj) 
with r^ = ©5 {v-'^) , rj = O5 (v-^) ■ 

Proof. Recall the definition of given by ([s]) and introduce the notation {t',v') = 
f{t,v). Define the two functions (see e.g. \28^ for a motivation of the formula 
defining z9) 

(8) d{t)^J'^J^r^{s)ds, I{t,v)^J r^{s)ds 

It is immediate to observe that i9(t') — ^{t) + I^^{t,v). Since the expression 
defining / is implicit, we find it convenient to use a suitable approximation in our 
computations. Define J : A \ 5 — > M as 

(9) 2j-'Ji;v)^ve + ei + l—. 

S V 

We claim that h : {t, v) (i9(t), J{t, v)) is the required change of coordinates. The 
first step is to obtain an approximate solution of ([2]). Since i is Lipshitz continuous, 
we can find the solution by iteration. Let (5t(o) = and define for n > 

eit) + iit + st^„_,)it,v)) 

6t(n)[t,V) = . 

V 

Then — <5i(„_i)|| = O (i'^") and thus — > St uniformly. Consequently, if 

we express the solution as 

(10) (5i = £ <5t„ with St^it, v)^^, 

n—l 

we can then find the functions a„ by the previous argument. In particular, outside 
R we obtain that 

(11) St,{;v)^-, St2{;v)^^e, St^i' , v) ^ ^ + U) . 

Assume now that (t, v) ^ R. By expanding Q in Taylor series and using equations 
(11 1 it is immediate to check that 

(12) J^I + 05{v-^). 
Recall that 

r^{t,v) - ^{t')-m - J{t,v)-\ rj{t,v) = J{t',v') - J{t,v). 

IS 



3.2 



Thus estimate ( 12 ) immediately yields r^ = ©5 (w ) . The proof of Lemma 
thus complete once we prove that 

J{t',v')~J{t,v) = 05 (v-^). 

We begin by introducing a convenient notation. Fix {t, v). Recall that rj = Jo/— J; 
denote J — J{t, v), J' = J{t', v'), i = £{t), £' = £{t') and likewise for all derivatives. 
Notice that Jv^^ is a polynomial in v^^ with coefficients given by smooth functions 
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of t. Using ( 10 1 we can express 6t in similar form, thus, by expanding in Taylor 
series the smooth function £ and its derivatives we can write 

where r*j = O5 {v^^) ■ It amounts to a simple but tedious computation to show 
that our choice of J implies 60 = 0, 61 = and 62 = 0. Here we will only sketch 
the main steps of the computation. First, we obtain an expression for Sv 

dv = v' - V = -2i' = 6vq + Svi + Sv2 + (v^^) where 



6vo 



-2£ 



Svi 



SV2 



Next, we expand rj in Taylor series and collect terms of order v or higher in the 



function r j 

rj = dtJ{5ti + 6t2 + Sts) + dyJ{5vo + 5vi + 5v2)+ 

+ ^9ttJ (Stl + 2StiSt2) + dtvJ{5tiSvi + 5ti5v2 + 5t2<5wi)- 
+ ^dtttJStl + ]^dttvJStl5vi + r*j. 
Using the explicit form (|9]) it is then simple to obtain 
5o = (.vSti + £5vq 

bi = ev^St2 + + U)vSti + ^iv^Stj + tvSvi + ivShSvo. 



and finally 



62 = iv^6t3 + + (.'i)v^5t2 + -{2e££ + f"e)v5h + 
+ iv^StiSt2 + ]^{U£ + ££)v^5tl + ^£v^5tl+ 

+ iv'^5v2 - \££5vq + Iv^ShSvi + £v^St2Svo + l^iv^StlSvo- 

Now it is possible to conclude by substituting bj into the formulae obtained previ- 
ously. □ 

Proof of Theorem [7} It is simple to check by inspection that (r, /) with t — Id are 
smooth coordinates on R for sufficiently large /. Let {t,v) G R, {t,v) — F{t,v) 
and {i,v) = f~-^{t,v). We use the convenient shorthand notation J = J{t,v), 
J = J(t,v) and J = J(t,v) and similarly for z9, and By iteration of Lemma 
13.21 we obtain 

J{i,v)-J{t,v)^05 {v'^). 

We then claim that 
(13) J-J^Ij{£+^£- 



£+ 

tv{l~—)-£o 

V 



1 ^£+ - £- 

+ T'^ — 

4 V 



{tv ~ £of - -£\ 
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In fact notice that 

i{t) =i+ + (v-^) 

i{t) = i+ + in + o {v-^) 



= r + ti+o (ii^^) 



i{t) = 4 + tn+ ^e+P + o i{t) = 4 + tt + ]^tP + o 



and moreover 

(14) f=t--^ + ^t-~^ 

V V V 

By the definition of J we thus obtain 



Oiv 



- 2£{t). 



{l(t) - l{i))v+ 

1 t{t)'i{t) ~ e{i)£{i) 



from which ( 13 ) follows by a straightforward computation. Notice that by definition 
we have 

2J-^J = v£o + i+vi+io£+ + O {f}-'^) 
2J-^J = Mo + £-vi+£ot + O . 

Next, by definition of 

t = J£l^{l + J£o£+d) + O {v-^) 

from which we obtain 

(tv{l + £+/v) - ^o) = 24(r - 1/2) + O (f-^) . 

Therefore we can rewrite J as follows 

J - J + A(f - 1/2) -^^{{t- 1/2)2 _ ^/^2) + O5 {J-^) . 



Using estimate (12 1 we thus conclude that 

I{lv)~I{t,v)^I{i,v) + 0^ 

We now prove that 

(15) r = r — J mod 1. 



By Lemma 3.2 and definition we have 

J§ = T-J-l 



On the other hand, using the definition of ?9 and the approximate expressions for 
J given above, we obtain 



2io 



From the last equation we obtain, using ( 14 ) , that 



Now (15) follows from (12). 



□ 
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4. Hyperbolic case. Properties of the limiting map 

The proof of ergodicity of the map F has been first given in |26J. Stronger 
statistical properties claimed by Theorem |2] follow from the following general result. 
Let G be a piccewise linear hyperbolic automorphism of and denote by S+ and 
S- the discontinuity curves of G^^ and G, respectively; let S = S-U S+. For any 
positive n e N let Sn = G"^^5+ and 5_„ = G^("^^^S'_; assume for convenience 
5o = 0; let 5(")=UL-„^fc- 
Proposition 4.1 (Chernov, [10^). Assume: 

(a) Si n Sj is a finite set of isolated points if i ^ j; 

(b) S is everywhere transversal to the invariant stable and unstable directions; 

(c) for every m > 1, the number of components of S^^^^ meeting at a single point is 
bounded by Km for some constant K; 

then G is ergodic, mixing and enjoys exponential decay of correlations for Holder 
observables. 

Proof of Theorem^ If Tr(yl) > 2 then F is a piecewise linear hyperbolic automor- 
phism of T^; we recall the explicit formula: 

F : {t,I) ^ {t - I mod 1,/ + A((r-/ mod 1) - 1/2)). 

Thus it is easy to check that 5_ is given by the diagonal circle r = / and S+ is 
given by the vertical circle r = 0. It is then a simple linear algebra computation to 
prove that the stable and unstable slopes are given by the solution of the quadratic 
equation — Ah + A = 0; thus we immediately obtain item (b) in the hypotheses 



of Proposition 4.1 Since dF is constant at any point where it can be defined, the 
n-th image of any line segment is a finite disjoint union of line segments parallel to 
each other; hence each point p G can meet at most two of such segments, which 
proves item (c). Finally, unless the initial line segment is aligned to an invariant 
direction (stable or unstable), the slopes of line segments belonging to images at 
different times are different, which proves item (a) and concludes the proof. □ 



5. Elliptic case. Growth of energy. 

Proof of Theorem In order to prove item (a) we will prove that for each v suf- 
ficiently large there exists a stable periodic point {t^,v^) G R whose orbit satisfies 
condition ([6]); stability of the fixed point then implies that ^ holds for any initial 
condition (to, vo) in a small ball around (<*, w*). We already noticed that the point 
(f, J) is a stable fixed point of the map if f = 1/2 and J S Z; in order to prove 
existence of a stable fixed point of the first return map F we would need to prove 
that the fixed point of F satisfies the non-degenerate twist condition. However, 
since F is piecewise linear, we actually need to consider the first return map as a 
O {J~^) perturbation of the map F = F + Fi and check that F satisfies the twist 
condition. Since the perturbation term is small up to derivatives of sufficiently high 
order, we can conclude. 

Fix once and for all J G N such that \J{0,v) — j| < 2; let A — exp(i0) be the 
multiplier at the fixed point (f, J) of F. Since A e (0, 4) we have A 7^ 1; then F will 
have a fixed point close to (f, J) that we denote by (f , J); introduce the parameter 
e — Ai//; by inspection is is easy to see J — J and f — t + e/(12A) + O (e^). 
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Introduce coordinates (cr, J) in a neighborhood of the fixed point (r, J) such that 
(r, J) = (f + (T, j + J). The expression for dF in these new coordinates is: 



dF{a, J) 



1 -1 
K + ea 1 — K — ea 



where k = A + e^/(6A); denote by A = exp(i6') the multipher of the map F at 
(f, J): it is immediate to check that 

(16) cos e^cos9~ eV(12/A). 

In order to check the twist condition we perform a complex change of variables 
(ct, J) H' (z, z) such that the map can be expressed as follows 

F : z H> z + Az + A^z^ + A^zz + Agz^. 

Then (see e.g. p^) we need to ensure that: 

Notice that from the fact that F is symplectic we obtain that — \A^\^\ thus 

there are two possibilities for the twist condition to fail; either ^3 = = or A 
solves the equation 

It is easy to check that the above condition is given by either = or cos d = —1/4; 
both possibilities can be prevented by avoiding special values of e, according to ( 16 1, 
possibly by choosing a different J. Therefore, we just need to check that A3 7^ 0: 
from elementary linear algebra we find that: 

z = cr+(l-A)J z = cr + (l-A)J; 

changing variables we obtain: 

A-1 , 1-A-A , A-1 

Which implies IA3I 7^ and concludes the proof of item (a). 

The proof of item (b) is analogous to the proof of the corresponding result 
obtained in [15], Section 3, and will therefore be omitted. □ 



6. Hyperbolic case. Measure of accelerating orbits. 

Note that part (a) of Theorem |4] follows from part (b), however since the proof of 
part (b) is rather involved we give a direct proof in this section. We expect Lemma 
|6.1| below to be useful for a wide range of mechanical systems. In particular. 
Theorem |4]( a) is a direct consequence of Theorem [2] and Lemma [6. 1| 

Let X be a Borel space and F be a subset of X x N containing {(a;, to) : to > m\ 
for some to. Let $ : y — >■ y be the map 

$(2:, to) — (0(a;, to), to + 7(0;, to)). 

Assume that $ is asymptotically periodic in the following sense. Denote Tk{x, to) = 
{x,m + k) and consider "ifi; — Tf^^^Tj.. Assume that there exist a map tp : X ^ X 
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preserving a probability measure /i, and a function 7 : X — > Z such that for each 
M for each function h supported on X x [0,Af] and each I we have 

(17) \\ho v|/[. ~ho \E''||i2(^) where ^{x,m) = {;i]j{x),m + ^{x)) 

and is a product of /i and a counting measure on Z. Denote (a;„, to„) = $"(xo, mo) 
and let 

E = {(xo, mo) : m„ -> +00}. 

Lemma 6.1. Assume that 

(i) tp is ergodic with respect to fj,; 
(^^) /x 7(a;)dM(a;) = _ 

(Hi) <I> preserves a measure v with bounded density with respect to /z; 
(iv) \\j{x, m)||L~ < K. 
Then V{E) = 0. 

Proof. By \i\ we know that conditions (i) and (ii) imply that Xy.l^Xy.1 
is conservative. That is for each subset Y of finite Jl measure the Poincarc map 
^ : y F is defined almost everywhere. Let Y = X x [{),K + I] where K is 
the constant from condition (iv). By Rohlin Lemma applied to ^ for each e there 
exists a set D,^ dY and a number such that 

'jlip.e) < £ and Jl{{x,m) : #'(2;, m) <^ ioi < I < L^) < e. 



In view of (171 and condition (iii) there exists a constant C (independent of e) and 
a number k{e) such that for k > fc(e) we have ^{dk^^) < Ce where 

flk,e = {{x, m) e TkY : ^ T^f^e for I G N}. 

Let 

= U {^k{l/n^)^l/n^[J^k(l/n^),(l/n^)j ■ 

Note that ^{^l) < 00. On the other hand if (x,m) e E then, due to condition 
(iv), its orbit Orb(a;,m) visits Ti-Y for all k except for finitely many k. Hence 
Orb(x,m) n (fi n -B) ^ 0. Accordingly it suffices to show that D E) = 0. 
However, by the foregoing discussion, the first return map $ : (il n i?) — > (fi n i?) 
is defined almost everywhere and by Poincare recurrence theorem for almost all 
(x,m) G ^ n E we have Orb(x,m) D X x {m} 7^ 0. Thus for almost all points in 
QnE we have {x, m) ^ E. Therefore n i?) = as claimed. □ 

7. Hyperbolic case. Time of deceleration. 

7.1. Plan of the proof. Here we prove Theorem |4]^b). The argument of this 
section has many similarities with the arguments in [5J 1141 116) so we just indicate 
the key steps. 

The proof relies on the notion of standard pair. A standard pair is a pair i — (7, p) 
there 7 is a curve such that I7I < 1 where I7I denotes the length of 7, 7' belongs 
to an unstable cone, |7"| < Ki, and p is a probability density on 7 satisfying 
1 1 In p| lci(^) < K2. We let \i\ denote the length of 7. We denote by the expectation 
with respect to the standard pair 



E^(yl) = / A{x)p{x)dx 

J-y 

and by the associated probability measure, that is, Pe{^) ~ ^ei^n)- 
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An easy computation shows that if / is sufficiently large on 7 then the standard 
pairs are invariant by dynamics, that is 

j 

where Cj = 1 and ij are standard pairs. We need to know that most of "fj in 
this decomposition are long. To this end let r„(a;) be the distance from Xn to the 
boundary of the component "fj containing Xn- 

Lemma 7.1 (Growth lemma). 

(a) There exist constants C > and < I such that 

P,(r„(x) <e)<Cs + Vf{ro{x) < er). 

(b) There exists a constant Sq such that if Uq > K\ In |7|| then 

Pi{rn{x) < £0 for n^no,...,no + k) < CO'' . 

The Growth Lemma is the key element of proving exponential mixing for F (see 
[TT1[S]) and the argument used to prove the Growth Lemma for F shows that this 
property is also valid for small perturbations of F . 

Given a point x let Ta be the first time /„ < a if /q > a and be the first time 
/„ > a if /q < a. Let Ta b be the first time then either /„ < a or /„ > b. 



The proof of part (b) of Theorem 
is an extended version of Theorem [5 
second one gives an a priori bounded needed to handle small velocities 

Fix < a < 1 < 6. Denote 



4] depends on two propositions. The first one 
It will allow to handle large velocities. The 



CO 



(18) 1)2^ y // A{x)AiF''x)dx 



n— — 00 



where F = G o Ta and G and Ta are defined by ([5| . 

Proposition 7.2. Let x be distributed according to a standard pair £ such that 
/ - /o oni and |£| > Then 

(a) The process 

converges to the Brownian Motion with zero mean and variance D^t which is 
started from 1 and is stopped when it reaches either a or 5; 

(b) There exists S > such that \E£{It^j^^ j_^^) — /o| < C'/q^''; 

(c) There exists 9 <1 such that Vi{TaioMo > kl^) < ma,x{0'' , 1^^"°); 

(d) LetT*f, = min[Taio,bio,lf))- Then 

¥i{rT'{x)<e)<Clle. 



Proposition 7.3. Given e > there exists K{e) > such that if \i\ > Iq th 

P,(Tc > Kie)Q < s. 
Note that Proposition 1 7 . 2| implies that 

ViiTsio > til) ^ > t) 



en 
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where Tf denotes the first time the Brownian Motion from Proposition 
S. Indeed 



7.2 



reaches 



hiTsi, > til) - Vt{Tsio > til and Tsi, < Tai„)\ < ViiTgi, > Tai„). 



By Proposition |7.2| the RHS can be made as small as we wish by taking A large. 
On the other hand by Proposition |7.2| 

MTsio > til and Tsi„ < Tai,) ^ V,{T^ > t and Tf < Tf) 

and the last expression can be made as close to P(r5^ > t) as we wish by taking A 
large. 

Next note that it is enough to prove Theorems Qb) and [5] with v replaced by /. 
Indeed, in view of ([12]), ^ and ([3]), we have 



(19) 



which shows that v can be replaced by / in Theorem [sj Also ( 19 1 allows us to 
squeeze the first time v goes below C between the time / goes below Ci and the 
time / goes below C2 and, in view of Proposition |7.3[ the times to go below Ci and 
C2 satisfy the same estimates. 

We are now ready to derive part (b) of Theorem [4] from Propositions 7.2 and 

ESI 



We have 



i{Tc < til) < fiiTsio < til) ^ IP(rf < t) 



and the last expression can be made as close to V{T^ < t) as we wish by taking 6 
small. 

Conversely 

MTc > til) < P,{Tsi„ >{t- K{e)S')ll) + ¥,{Tc{xtsJ > K{e)S'll) 



where K{e) is given by Proposition 7.3 The first term can be made as close to 
P{Tif > t) as we wish by taking small S. To estimate the second term note that 

PdTc{xTsJ>K{e)Sll) 
= P(rc(xT„J > Kie)S'll and rT„„(a;) < {SIo)-'°°) 
+ViTcixT,J > K{e)5^ll and rT,,^{x) > (Wo)"'™) = / + J- 

Next 

I<F{rrs,,{x)<{SIo)-'"''))^0{I^'') 
by Proposition |7.2[ d) and 

II^P{rT,,„{x) > /o-i™))P,(Tc(xT„J > K{s)5'll\rT,,^ > (SIo)-'"')) 

< ¥i{Tc{xTsJ > Kie)6^ll\rT,,^ > (SIo)-''')) < e 

where the last inequality follows by definition of K{e). 

This completes the derivation of Theorem |4](b) from Propositions 7.2 an d|7.3[ It 
remains to establish the propositions. Proposition |7.2| is proven in section |7.2| and 
Proposition |7.3| is proven in section [73} 
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7.2. Central Limit Theorem. Let 

Fj{I, r) = F{I, t) + [/]-i(0, Ai((f - l/2f - 1/12)). 

Note that f\ approximates F up to error 0(/^^). Next consider a mapping of the 
given by 

F^(/,r) = i?(/,r) + iV-i(0,Ai((f- 1/2)2 _ ^/i2)). 

Then F\ locally covers Fj/j; also Fjv preserves the measure d/dr. The proof of 
Theorem [2] shows that Fjq is exponentially mixing. In particular, if \l\ > Eq then 



¥.i,{AoF^)^ Ad/dT + 0(6i"). 
We use this property to establish the following estimate 

Lemma 7.4 (Averaging Lemma). Suppose that |£| > Let n ~ KIwIq where 

K is sufficiently large. Let A he a piecewise smooth periodic function. 

(a) EeiA o F») = //^, AdLdT + 0(V'+*); 

(b) There is L > such that 



Ee{A{F''x)A{F"+''x)) = jj A{x)A{F''x)<:\IdT + 0{i;^^L^) 

The proof of this lemma is similar to the proof of Proposition 3.3 in [7]. The 
proof of part (a) proceeds in two steps. First, if \£\ > Eq then we use the shadowing 
argument to show that 

(20) EeiA o F") = E,{A o F^'J^,) + 0{Lo''+') 

and then use exponential mixing of F^ig]- In the general case we find a function 
n{x) < ^ In/o such that 



E,(A(F"(-)x))=^c,E,^.(A) 



and J2\ej\<eo '^i — ^0 ^'^"^ then apply (20) to all long components £j. 
To prove part (b) we first use the foregoing argument to show that 

E^(yl(F"a;)A(F"+'=a;)) = // A{x)A{F''x)dLdT + 0{L'^-^L'') 



(the factor L'^ accounts for the exponential growth of the Lipshitz norm of A{A ' 
F'^)) and then use the shadowing argument again to show that 



A{x)A{F''x)dLdT ^ // Aix)A{F''x)dLdT + 0{L^ 



It is shown in Appendix A that Lemmas 7.1 and 7.4 imply parts (a) and (b) of 
Proposition |7.2| We note that the error bound 

0(/-(2-<5)) ^ j-i 

is needed to compute the drift of the limiting process; to compute its variance it is 
enough that F = F + 0{L~^) and that F covers F which satisfies the CLT in the 
sense that 

Normal(0, D^) 
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where the diffusion coefficient is given by tlie Green-Kubo formuia (18 1. (In 
fact (18) is the Green-Kubo formula for F = GFG^^ but F and F clearly have the 



same transport coefficients.) 

Next, part (a) of Proposition 7.2 implies part (c) with k — 1, that is, there is 
9 < 1 such that 

(21) MTaloMo < lo) < e. 

For fc > 1 we argue by induction applying (21 ) to all long components of F^'^~^'>^oj 



which have not escaped by the time (fc — 1)1^ ■ Finally 

ViiKb <^)< E ^drmix) < e) 

m—K In Iq 

SO part (d) follows from part (a) of Lemma |7.1[ 

7.3. A priori bounds for the return time. Let (Tq be the first time when — 
2™"! < A. For j > 1 we define aj inductively as follows. Assume that crj-i was 
already defined so that |/o- _i — 2"'^^^ | < A. Let aj be the first time after aj-i when 
either |/„ - 2'"^+i| < A or \I„ - 2"^-i| < A. Let cr, = min(o-i, cr,_i + 23"^-i). 



-lOOm 



j 2™3 < / then we stop otherwise 



If either > aj^i + 23™j-i or r^. < 2" 
we continue and proceed to define crj^i. If we stop we let j* = j be the stopping 
moment. If we stop for the first or the second reason we say that we have an 
emergency stop, otherwise we have a normal stop. By the discussion at the end 
of section [TTT] the lower cutoff in Proposition |7.3| is not important so to prove the 
proposition it is enough to control the first time when /„ is close to 2™ with 2™ < /. 
In other words we need to control aj*, especially if it is a normal stop. Also since 



(To is unlikely to be large by part (c) of Proposition 7.2 (in fact, part (a) would also 
be sufficient for our purposes) we need to control CTj* — ctq. 

Let ¥j be the cr-algebra generated by (mo,cro)-''' : ('^^j ? "'j ) • Proposition 7.2 
implies that 



1|F, 



P^(m,+i = 
j+i is an emergency stop|Fj) = O (2"9'^™3^ 



- 1|F,) 



1 

3 
2 

3 
O 



o(l), 



0(1), 

-97m, 



Let be a random walk with — f^o and 

P(^,+i = + 1) = 0.4, 
Let A^- be iid random variables independent of ^s such that 



P(0+i - - 1) = 0.6 



V{Aj = fc) 




if fc > fco 
otherwise 



where fco is sufficiently large and K — jjt^. Let Aj = min(2^^J Aj, 2^^^). Proposition 
7.2 allows us to construct a coupling such that for j < j* 



1 ' 



< O-o + E ^r. 



m—0 
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Now a standard computation with random walks shows that Proposition |7.3| is valid 
for the random walk itself. Consequently, given e, there exists K{e) such that 

Unfortunately j* need not to be a normal stop, it can be an emergency stop as 
well. To deal with this problem let 

Pfc — is an emergency stop and mj- ~ k). 

Denote ^ki = {j* is an emergency stop, vrij* = k and j* is the Z-th visit to fc}, 
Vki — {k is visited at least / times }. Then 



Pfc < ^P^(r!fcO = ^PfC(4^)IP^(^^H|^H)■ 
By Proposition |7.2[ d) 



1=1 1=1 



while the existence of the coupling with the random walk discussed above implies 
that 

Vi{Vu)<0'. 

Therefore 

(22) pfc < C2-9^^ 

Accordingly, by choosing / large enough we can make the probability of an emer- 
gency stop less than 0.1. However we can not decrease that probability below e/2 
if I is fixed, so more work is needed. 

First, we note that (wj- > mo/2) = 0(/q"^°) so it can be neglected. Secondly, 
an argument similar to one leading to ( 22 ) shows that 

P.K. </o-^"") = o(v^°). 

Next, if r„., > nij* < mo/2 and j* is an emergency stop let a be the first 

time after ctj. such that rg{x) > £q. By the Growth Lemma 



e{a - aj* > KIhIq) < Iq 



100 



if K is large enough. 

If (T — aj* < Klnlo then we can repeat the procedure described above with x 
replaced by x^- If the second stop is a normal one we are done, otherwise we try 
the third time and so on. We have 

P£(First k stops are emergency stops) < (0.1)'° 

which can be made less than e/10 if k is large enough. Next we have 

Vdk* <k,Tj<Kie)I^)<^- 

5 

since the first try takes less than ^^-^Iq with probability greater than 1 — ^ and 
all other tries take time 0{Iq) since with overwhelming probabili ty w e start those 
tries below level 0{\/Tq). This concludes the proof of Proposition 7.3 
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8. Hyperbolic case. Dimension of accelerating orbits. 

Proof of Theorem Foliate the phase space by hne segments paraUel to the un- 
stable direction of the limiting map F. It sufBces to show that, given s < 1 there 
exists / such that if F is a leaf of our foliation and / > / on 7 then HD(F OE) > s. 

By Theorem [2] the limiting map satisfies CLT. That is, for any unstable curve 
7, if the initial conditions are distributed uniformly on 7 then converges to 

a normal distribution with zero mean and some variance D (here we are using the 
notation F^^tq^Iq) — (f„,/„)). In particular there exists a constant k > such 
that, for sufficiently large no, we have 

(23) Pf (/„„ -Io> > ^ 

where £ denotes the standard pair (7, Const). Moreover, given 5 < (5, we can find 



no so that (23) holds uniformly for all curves of length between S and 6. Let f„(a:) 



denote the distance from F"-x to the bound ary of the component of i^"F contain- 



ing F"x. By the Growth Lemma (Lemma 7.1) if S is sufhcicntly small than for 
sufficiently large no we have 

P,(f„„(x) <3<5) < ^ 

provided that 7 is longer than S. By Theorem [l] we can take / so large that if /o > / 
on 7 then 

Pf(/„o - Iq > Ky/n^ and r„o(x) > 3(5) > - 



where F{Iq,tq) — (/„,r„) and rn{x), as before, denotes the distance from F"x to 
the boundary of the component of F"F containing F"a;. Note that any curve of 
length greater than 36 can be decomposed as a disjoint union of curves with lengths 
between 6 and 25. Hence F"°'y D IJ^ where on each 7^ the action grew up by at 
least kuq and the total measure of IJ^- F^^'-°jj is at least mes(7)/5. Next, suppose 
that (5 < I7I < 26. Then we have |i^^"°7j | > 2{\+e)"o l^l number of curves 

is at least j^^, (A — e)"" where A is the expansion coefhcient of F and e can be made 
as small as we wish by taking / large. 

Continuing this procedure inductively we construct a Cantor set inside F such 
that each interval has at least jq{^ — e)"° children and ratios of the lengths of 
children to the length of the parent are at least 2{\+e)^o ■ follows that the resulting 
Cantor set has dimension at least 

In (3^, (A -£)"») 
In (2(A + £)"«) ■ 

This number can be made as close to 1 as we wish by taking no large and then 
taking / large to make e as small as needed. □ 

Remark 8.1. The Cantor set above is constructed by taking as children the sub- 
interval where energy grows by n^/n^. However, the same estimate remains valid 
if we take sometimes children with increasing energy and sometimes children with 
decreasing energy as long as /„ always stays above /. For example we can require 
that the energy grows until it reaches 21 then decays until it falls below ^ then 
grows above 3/ then decays below ^ then grows above 4/ etc. Then the argument 
presented above shows that the set of oscillatory orbits has full Hausdorff dimension. 
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We expect that this set has also positive measure but the proof of this fact seems 
out reach at the moment. 

9. Conclusions. 

In this paper we considered piecewise smooth Fcrmi-Ulam ping pong systems. 
Near infinity this system can be represented as a smaU perturbation of the identity 
map. Small smooth perturbations of the identity were studied in the context of inner 
pi] and outer (HZj) billiards. In this case, after a suitable change of coordinates, 
the problem can be reduced to the study of small perturbations of the map 

Tn+l — Tn + Uj{Jn), J71+I = Jn- 

This map is integrable so the above mentioned problems fall in the context of 
small smooth perturbations of integrable systems (i.e. KAM theory). In the case 
of piecewise smooth perturbations the normal form also exists: it is a piecewise 
linear map of a torus. However in contrast with the smooth case the dynamics 
of the limiting map is much more complicated and, in fact, it is not completely 
understood, especially then the linear part is not hyperbolic. In this paper we 
described for a simple model example: 

(i) how to obtain the limiting map and 

(ii) how the properties of the limiting map can be translated to results about the 
diffusion for the actual systems. 

However, there are plenty of open question on both stages of this procedure. For 
example, for piecewise smooth Fermi-Ulam ping pongs it is unknown if there is a 
positive measure set of oscillatory orbits such that 

lim inf u„ < cx) , lim sup w„ = oo 

in fact no such orbit is known for A € (0,4). This demonstrates that more effort 
is needed in order to develop a general theory of piecewise smooth near integrable 
systems. 
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